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ABSTRACT 


‘An investigation of using the response of a €1nclvt at. selected test 
frequencies to isolate faulty circuit components ts made.. A procedure 
using a sensitivity approach for -intelligent selection of test frequencies 
is developed. The developed procedure is tested and the results are 
compared with results using conventional procedures.. The effect of 
random, within tolerance variations of nonfaulty components on the 


results is studied for both conventional and developed procedures. 
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I. INTRODUCTION 


Maintenance of electronic equipment is a problem that increases as 
equipment becomes more complex. With ever increasing costs (of equip- 
ment) it also becomes imperative to decrease equipment downtime as 


much as possible. To this end a fast, reliable method of locating 


. circuit failures must be used. 


“It has been customary to isolate failures by testing individual compo- 
nents. This is usually done by isolating the fault to a particular stage, 
followed by subassembly isolation and finally component isolation by 
individual testing. 

With the availability and flexibility of digital computers, a pos- 
sible solution to the problem of excessive downtime is to use a computer 
to help isolate the fault. The first step in accomplishing this is the 
development of a testing and isolation procedure which could be progran- 
med for the digital computer. 

A recently presented procedure has em araianed by a number of 
individuals [1], [2], [3], and [4]. The procedure was described by Seshu 
~ and Waxman [1] with other investigators testing and making applications of 
the procedure with few modifications. 

The basic procedure is to examine the response of a circuit under test 
at selected frequencies. Comparisons of the results with a list of “library 
results" can pinpoint a faulty component. 

Most investigators have used test frequencies near the poles and zeros 
of the function under consideration. Other frequencies have been chosen 
between the pole and zero frequencies [2] and [4]. In this study an 
alternate method of test frequency selection was developed and tested and 


results were compared with results using other methods. 
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The procedures for signature comparison used by other investigators 
have not been explained in detail in most reports. However, it appears 
that a nearest neighbor approach is the most common technique used to 
date. Two alternate methods were examined and are discussed with 
conclusions as to their practicality. 

Previous methods of fault determination after signature comparison 
have also not been documented. Two possible methods were tried. The 
first method used computer selection for each network function and manval 
correlation for fault isolation. The second method. was completely program- 
med for digital computation. | 

As far as can be determined from available reports no investigator 
approached the problem of error due to in-tolerance variations of com- 
ponents. This study melded an investigation of that problem with both 


the conventional and the proposed procedures. 
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TI. PROCEDURAL BACKGROUND 


A. FUNDAMENTAL CONCEPT 

All components of a given network play a part in the formation of any 
of the network functions. The sensitivity of any given network function 
to a particular component is generally unique to that component. That is, 
in a network with K components a network function will have K sensitivity 
functions, one for each eoTpenenies A change in the value of a particular 
component may cause a decrease in the magnitude of a network function at 
one frequency and an increase in the magnitude at a different frequency. 
There may in fact be frequencies at which the function sensitivity is 
near or equal to zero for a particular component, excluding catastrophic 
failures. 

Specific knowledge about this behavicr can be used to detect and 
isolate component value changes. If, for example, it is known that the 
input impedance is insensitive at low frequencies to changes in a specific 
component, one would assume that component was not faulty if the input 


impedance exhibits a significant change at low test frequencies. 


Dee ePOLE—shl El TECHNIQUE 
The component variation effect on frequency response can be observed 
with the pole-shift technique. A network function can be represented by 
a constant and a unique set of poles and zeros; i.e., 
sed (Sez los2)) ee Sa2 7) 


where 


4 
o—~ 
¥) 
batt 
It 


any network function 


aS 
if 


constant 


1] 


Z. = ith zero frequency 


i 
-t 
j h 


in 


pole frequency. 


Every component plays a part in pole and zero locations. However,. a 
few components may dominate the location of a given pote or zero frequency. 
Thus, if one could detect the shift of the poles and zeros of a given 
function, it may be possible to determine component variations. A simple 
example can best illustrate this point. 


Example 1 - Consider the network shown in Figure 1. 





Figure 1. Example 1 - Low Pass Filter 


The network function 7 (s) is given by: 
] 


We feo een 


] 5° + 10s + ] 


The pole and zero shifts were determined by varying each component 
value by ten percent and observing the corresponding critical frequency 


shifts. The results: 
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Location (radians/sec) | | Percent Shift 


Bait ion Py P2 a @ Py ae 72 

Nomina] .101 9.899 1.113 9.887 0 9) 6] Q 
Ry variation  .10] 9.899 1.020 9.889 0 0 -62 02 
C, variation  .092 9.908 1.010 9.899 8.9 -Q9 Seas) Se 
L, variation .101 8.990 1.114 8.977 . Q 9.2 .08 9.2 
R, variation .092 10.908 1.101 10.899 8.9 10.2 R.O7e 1022 


TABLE II-1. Percent Pole and Zero Shifts 


It is immediately obvious that a change in Re or C, will cause a greater 
shift in Z, than will a corresponding change in L, or R.- Also, Po and Z5 


are most affected by L and R, and Py is affected only by C, and Ro. 


C. CHANGES IN RESPONSE FROM POLE SHIFT 

In practice the exact location of the critical frequencies can usually 
not be determined directly from the circuit. However, this does not rule 
out the pole shift technique. Any shift in critical frequency causes a 
change in the magnitude of the given function at all frequencies greater 
than the shifted one. This can be seen from the simple Bode representation 


of a pole: 


1G (F) | 





Figure 2. Bode Plot of Simple Pole 
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If os is changed to le (indicated by dashed line), the function 
magnitude, at any frequency, fo. greater than Fos is changed (increased 
in this case). The magnitude is not affected at frequencies less than 
hens 


C 
Now, from the network in Figure 1] we obtain: 


2 10 : 10 


ce gm ee Eee = 


1 s*4#1is +11 (s + 9.887)(s + 1.113) 


Figure 3a is the Bode plot for the function Vo/V) - Figure 3b and 
3c show the effect of increasing Ry and of increasing Lys respectively. 
The nominal response is shown as a dashed line. 


H( = .099 Py = 1.57 Hz 


Wap * 





01 5 i: 10. t 


Figure 3a. Vo/V) Bode Plot. 





1 ie f 


Figure 3b. Vo/Vy » Ry Increased. Figure 3c. Vo/V) » L, Increased. 





It can be seen that an increase in Ry Causes a decrease in the 
magnitude of the response between P, and Po: A lesser decrease is seen 
at frequencies greater than Po- The response between Py and Po 1s 
increased for an increase in Ly while the response above Po is decreased. 
Thus, if one examines the response at a frequency between Pp and Po and 
at a frequency greater than Po it will be a reasonably simple task to 
distinguish between an increase in R or L.. This example is an over- 
Simplification of the procedure, but the application of which can be 


appreciated. 
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TIT. GENERAL PROCEDURE 


A. BASIC PROCEDURE 

The basic procedures outlined by Seshu and Waxman [1] was followed 
with modifications as explained herein. Two network functions were used 
concurrently. Work done by Maenpaa, Stehman, and Stahl [5] indicated 
that the redundancy from using two functions was more than adequately 
compensated for by the increased fault detection capability. 

The symbolic network functions were calculated by hand. It should 
be pointed out that computer programs exist which calculate symbolic net- 
work functions. These programs should be used for more complex circuits. 

The functions were numerically evaluated and factored to determine 
the poles and zeros. Test frequencies were then chosen. Previous 
procedures generally selected test frequencies based on the pole and 
zero locations. As a test of the procedure several test frequencies were 
chosen between the critical frequencies as well as O hertz and one 
frequency above the highest critical Freauency: The actual selection 
was programmed for computer selection and is described in Section VII. 

It was felt that the procedure of selecting test frequencies based 
on poles and zeros was, although intuitively satisfying, a somewhat 
illogical procedure. The concern was not with the function itself but 
with the behavior of that function witn a parametric change. Thus, 
sensitivity functions were calculated for each element. Based on these 
functions, six test frequencies were chosen. This procedure is discussed 
in Section IV. 

After selection of test frequencies, a worst case response, for 


component values within 10 percent tolerance limits, was computed for 
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each test frequency. This defined a nominal range of performance. This 
is not necessarily the best method for definition of nominal range. This 
range might be specified for a particular-circuit. However, in the 
absence of such specifications the above procedure was followed. 

Each component was then allowed to take on several discrete values 
out of tolerance limits. All other components were held fixed at their 
nominal value. The response at each test frequency was then calculated 
-and catalogued for the entire collection of sets of component values. 

The catalogued response Ace were then quantized into nine levels 
and labeled one through nine. Quantization level five was reserved for 
response within nominal limits. The setting of the levels is explained 
in Part B. The quantized responses for a given set of component values 
is defined as the signature for the given circuit configuration. 

The signatures were computed for the following cases: 


value = zero 
+ 20% 
+ 50% 
100 times nominal. 


The entire set of signatures for a given function is defined as the 
library for that function. There is no theoretical limit. on the size 
of the library. However, the inclusion of many more circuit configurations 
would probably only complicate the matching procedure and not increase 
detection capability. 

Using the primary test circuit the input impedance and voltage gain 
at each test frequency was measured. Then, after quantization, the signa- 
tures measured were compared with the library signatures and a "closest 
match" was made. Thus, the faulty component was identified. 

In order to test the procedure, several error conditions were set and 


the response calculated. In order to more closely approximate actual 
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circuit conditions, all component values were allowed to vary randomly 


‘(uniform distribution) within 10 percent tolerance limits. Calculations 


were made with and without the random variations to determine how much 
effect this would have on the results. 


The results of all tests are presented in Section VI. 


B. SETTING OF QUANTIZATION LEVELS 

.With the exception of selection of test frequencies the most critical 
Step in the procedure 15 he setting of quantization levels. Two methods 
were used successfully. One was programmed as part of the main analysis 
program. The second method, although not programmed is possibly adaptable 
to computer selection. Both methods have nine levels of response with 
level five reserved for the nominal range. 


The nominal range was selected by taking the range between the 


maximum positive and maximum negative deviation from normal for any 


element value at + 10% from normal. The programmed method then set the 


other-eight levels as shown in Table III-1 in which F is a multiplying 
factor used to force the possible response range to cover all quantum 


levels with levels one and nine assigned only for extreme cases. 


Quantum Level Normalized Response 


1 (o, 1 + 13 FN] 

2 (1 + 13 FN, 1 + 7 FN] 
3 (1 + 7 FN, 1 + 4 FN] 
4 Csr Ne 1 +N] 

5 (1 +N, 1 - M) 

6 [l - M, 1 - 4°FM) 

7 fl - 4 FM, 1 - 7 FM) 
8 [1 - 7 FM, 1 - 13 FM) 


LO 


[1] - 13 FM, 0) 
TABLE III-1. Quantization levels. 
18 


The second method used a more inituitive approach. The maximum 
deviation levels, 1 and 9, were selected based on the maximum deviation 
from normal under any condition. The mid levels were then selected so 
that the catalogued or library signatures would cover all quantum levels 
as uniformly as possible. 

The manual approach is acceptable for small circuits but gets 
unmanageable quickly as circuit size increases. Parts of the procedure 
could be easily programmed. However, a suitable algorithm would have to 


be developed prior to complete programming. 


C. TWO COMPONENT VARIATIONS 

It cannot be assumed that only one component at a time will undergo 
Changes in value. ~In order to investigate the possibility of two compo- 
nent values changing simultaneously, two possibilities were considerec. 
If one component masks function changes due to variations in another 
component (i.e., no change in function magnitude), the entire procedure 
will not work. However, if one is primarily concerned with circuit 
performance, the failure to detect component variations is of no concern 
if there is no degradation in circuit performance. 

The second situation considered was one in which there was an obvious 
change in circuit response due to two components varying. The original 
procedure made allowances for varying only one component at a time. How- 
ever, a Slight modification allowed for the generation of a set of library 
signatures for two component variations. The set is presented in the 
computer output section. 

It soon became apparent that in order to use the two component 
variation library, one needed apriori knowledge of the number of faulty 


components. Since this is not usually available, the usefulness of two 


ie 





component signatures is doubtful. Additionally, an examination of the 
library signatures revealed that one of the two components varied usually 
dominates the response. Thus, if two components have failed the single 
component library would probably isolate the dominant component. After 
replacement of that component further testing could reveal the second 
faulty component. Consequently no further investigation of multi- 


component variation was undertaken. 
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IV. PRIMARY TEST CIRCUIT 


The primary circuit used throughout lhe investigation was a six— 
--element low-pass filter. (See Figure 4.) The unscaled frequency 
response, represented by the voltage ratio transfer function, Vo/V->: 
is shown in Figure 5 and the unscaled input impedance as a function 
of frequency is shown in Figure 6. On Figures 5 and 6 the poles are 


- -represented by X and the zeros by 0. 


Cy L, 


Figure 4. Primary Test Circuit. 


The component values were magnitude scaled by 100 and frequency 
‘scaled by 10°. This yielded component values of approximately the 
same magnitude. Ideally, the scaled cutoff frequency should be T Hz.. 
However, mien this scaling ie = .02 Hz and no difficulty was encounter- 
-ed. 


The scaled network functions of interest are: 
Vo 1634 


YW s'+17.675°+6.4625¢+2.601s+.2124 


e - 1634 
~ (st17.31)(s+.0997) (st1 .325+53.25) (st+] .325-j3.25) 


2] 








001 . 01 7 F 


Figure 5. Vo/Vy> Primary Test Circuit 


10 x se) Sx O 


001 01 5 Je 10. 


Figure 6. Vi/1)> Primary Test Circuit 
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My) 4417.675946.4625¢4+2. 6015+.2124 
mee es 22-6015t.2124 
I, -5.9825°4+.1475¢+.817s+.01634 





Sess. 0997) (s#1. S20 tI Ox So aaE B25 10425) 
2373 


The poles and zeros in Hz are listed in Table IV-1 for convenience. 


V,/1, V,/1, 
Zeros* (Hz) Poles (Hz) 
0518 (Second Order) = .0594 (Second Order) 
.0159 .00319 
Zeal 
* Zeros of Vi/T, are the poles of V/V). 


TABLE IV-1. Pole and Zero Magnitudes. 
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Vem oEeECiun Or Mest FREQUENCIES 


A. CONVENTIONAL METHOD 

The conventional procedure was followed with Beant to selection of 
test frequencies in order to test the method. It soon became apparent 
that the number of test frequencies was too large to handle with ease. 
Thus, only frequencies between the pole and zero frequencies were 
selected as well as 0 Hertz and one above the maximum critical frequency. 

The computer program written for implementation of the procedure (see 
Program 1) did the frequency selection automatically. The test fre- 


quencies were selected by taking a4] 


F 
f ite. “fy 
Lest Cie 2 
for every a where fe -: jth critical frequency. If 
at i Vaz i+] 
OF - f ; - f 
ft ne Ss ee ae ee E 


g C 2 
one of the frequencies was eliminated. Test frequencies were computed 
for the input impedance because the critical frequencies for Vo/V, are 
included in the critical frequencies for Vi/1,. One test frequency was 
eliminated due to the above inequality. This resulted in the selection 


of the following seven test frequencies for a fourth order filter: 


0.0 Radians 
001 Radians 
00498 Radians 
196 Radians 
5488 Radians 

8.49 Radians 

25.0 Radians 


TABLE V-1. Conventional Method Test Frequencies 
24 
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B. SENSITIVITY APPROACH TO FREQUENCY SELECTION 

Inspection of the signatures generated using the eonreneionell method 
showed that redundancy existed to the extent that the method was some- 
What inefficient. Two test frequencies, D.C. and the lowest nonzero test 
frequency, gave nearly identical results. This indicated one of the two 
frequencies was probably not required. It was also felt that unless 
further investigation was undertaken it would be possible to ignore a 
range of frequencies which could be extremely helpful. 

This investigation took the form of calculating and plotting. functions 
which defined the sensitivity of Vo/Vy and 1,/V, to variations in each 
component. These plots are presented as Figures 7 through 18. Figure 19 
1S a superposition of the plots of the sensitivity function for I,/V, and 
Figure 20 is the same for Vo/V,- Ideally VAY, Should have been used. 
However, the computer program used in this portion of the study did not 
allow for using Vi/1,- It was felt that no loss of information would 
result. 

The superposition plots Figures 19 and 20 are very revealing. The 
first observation is that Vo/Vy 1s not sensitive to changes in Ry and Ro 
at high frequencies nor to C's or L's at low frequencies. While this 
information is not new the sensitivity functions assist in determining 
when C's and L's have little effect. This can be used to eliminate the 
redundancy observed using the conventional method of frequency selection. 

The mid-frequency range is the most aneereccine range for the selection 
of test frequencies. For example, near .02 Hz Vo/Vy is very sensitive to 


changes in R and Cos less sensitive to changes in Las even less sensitive 


* The program NASAP-70 was used to calculate the various sensitivity 
functions. 
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to changes in Ro and C, and virtually insensitive to changes in L,. 
There is a peaking of the sensitivity functions for oa Ca and L., near 
.06 Hz. At this same frequency the sensitivity functions for Ra» Ro» 
and L, are near zero. 

The I/) and V, plots show that at near .3 Hz only the sensitivity 
function for Ry is high. The functions for L, and C, are low and the 
other functions are zero. At 3 Hz 1,/V, is about equally sensitive to 
changes in L. and Ry while the other functions are very close to zero. 
It can also be seen that the frequencies discussed in connection with 
Vo/V) are reasonably good choices for L/V). 

Based on the above observations the following six frequencies were 


selected as test frequencies: 


Number Frequency 


.00048 Hz 
102) ahz 
.06 Hz 
»299 Hz 
2.99 Hz 
49.97 Hz 


HO OF > Ww PY —_ 


TABLE V-2. Sensitivity Function Test Frequencies. 


After some preliminary work it was observed that the response at .06 
Hz was not well behaved. Since four sensitivity functions peaked for 
1,/V) random variation within tolerance limits of these components were 
being observed as significant changes in response. In order to overcome 
this problem, test frequency three was changed to .055 Hz. The new fre- 
quency was off the peak enough to be useable and still yield the desired 


information. 
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Figure 7. Sensitivity of Vo/¥, to RI. 
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Figure 8. Sensitivity of Vo/Vy COullia 


28 


1: 





4, [F.(w) | 


001 


SS 


01 si Ls 


Figure 9. Sensitivity of Vo/N, LOrGL. 
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001 


01 a 1; 
Figure 10. Sensitivity of Vo/Vy to R2. 
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Figure 11. Sensitivity of Vo/Vy to. EZ 
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.001 


01 a 1. 
Figure 12. Sensitivity of Vo/V, C6, C2. 
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.001 


-O1 .l 1. 
Figure 13. Sensitivity of 1 /V, to RI. 
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Figure 14. Sensitivity of 1,/V, ios! Gl 
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Figure 15. Sensitivity of 1/V, BOn Gly. 
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Figure 16. Sensitivity of I,/V, to Ro 
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Figure 17. 


a ie 
Sensitivity of 1, /V, COZ. 
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Figure 18. Sensitivity of I, /V, EQUI 2: 
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Figure 19. 
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1 /V, Composite Sensitivity Functions. 
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Figure 20. V/V) Composite Sensitivity Functions. 





C. COMPARISON OF METHODS 

The comparison of the results using the two methods of frequency 
selection is included in Section VI. Some comments as to applicability 
- should be made at this point. | 
If one computes sensitivity functions by using the derivative of 


the network function with respect to the element under consideration it 





becomes 
f(s) = de, 03) 
q(s) . ts) 
: ciel 
where: 


fo(s) = sensitivity function 


d/de, = derivative with respect to element k 


“— 
— 
tt 


ratio of polynomials representing the network function. 


Thus, the poles of the sensitivity function are the same as the poles of 
the original function. Consequently, if one selects test frequencies on 
the basis of function poles the selection is linked ta the sensitivity 
functions. Generally, however, the zeros of the sensitivity function are 
unique and are thus ignored by any selection procedure based on the poles 
and zeros of the original function. 

From the preceding paragraph, some of the advantages of using sensi- 
tivity functions for test frequency selection become apparent... However, 
that method has the disadvantage of requiring the calculation af sensi- 


tivity functions for each element. 


4) 





With the small circuit used as an example the value of sensitivity 
functions are demonstrated. With more complex circuits it appears more 
probable that important test frequencies might be overlooked if the sen- 


sitivity functions are not used. 
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VI. RESULTS WITH PRIMARY CIRCUIT 


A. RESULTS USING BASIC PROCEDURE 

Using the procedures outlined in Section III-A and V-A a set of 
signatures was generated for the six selected circuit conditions. 
These signatures are presented in Figure 25. A second set of signatures 
was then generated for the same circuit conditions with random variaticn 
allowed for the circuit elements. The second set of signatures was input 
to the selection program to determine what errors resulted from the random 
variations. 

Table VI-1 lists the input condition and the program output. The 


numerical result was determined as follows: 


Selection was correct condition 

Selection was correct element 

One of multiple selections was the correct condition 
One of multiple selections was correct element 
Incorrect selection. 


Ot & GW NH —_— 
i 


By correlating the results from VA/V, and Vj/1, aceite, MlSSss.0K 
“possible hit" rating was assigned for each test Signature as follows. 
A hit was assigned if both Vo/Vy and Vi/1, yielded a 1 or 2 or if one 
was 1 or 2 and the other was 3 or 4. A 3-5, 4-5, or 5-5 combination was 
called a miss. All other combinations were.called possible hits. For 
some cases deviations were warranted. For example, a 3-4 selection could 
have been called a miss if a component other than the given component 
was clearly indicated as being faulty. 

The result was 21 hits, 9 possible hits and 6 misses for 58% hits and 
83% hits or possible hits. Since only the circuit conditions used for 


Signature generation were used to generate the test signatures 1002 


43 





Circuit Program Numerical Program Numerical 
Condition Condition Result Condition Result 
Rl: Short R1: Short ] Rl: Short 1 
-50% R1: -50% 1 L2: Open 5 
-~20% Ri: Open 2 L2: +50% 5 
+207 L2: +20% 5 R1: +202 ] 
+50% R1: +20% 2 R1: +20% 2 
Open Rl: Open q L2: Short 
Rl: Open 3 
Li: Short L1: Short ] L1: -50% 2 
~50% L1: -50% 1 L1: -50% ] 
-20% L1: -20% 3 Ll: -20% ] 
R2: +20% 
+20% L1: +20% ] Ll: -20% V4 
+50% L1: +20% 2 L1: -20% 2 
-~Open L1: Open ] R27 4207 5 
Cl: Open Cl: Open ] Plc toes 5 
-50% Cl: +20% 4 Cl: -50% ] 
C2: -50% 
-20% C2: +20% 5 R1: -20% 5 
+20% l= 207, 5 L1: -50% 5 
C2: +20% 
+50% CZ e207, 5 Cl: +20% 3 
C2: +20% 
Short Cie Short ] Cl eShore ] 
R2: Short R2: ~-50% 2 R2: Short ] 
~50% R2: -50% ] R2: -50% 2 
R2: -20% 
~20% R2: -20% 3 Rl: -20% 5 
L2: +20% 
C2: +20% 
+20% G2: +20% 5 Cl: +20% 5 
C2: +20% 
+50% R2: +20% 2 R2: +20% 2 
Open R2: Open 1 C1: Open 5 
TABLE VI-1. Basic Procedure Results. 
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TABLE VI-1 (Continued) 


>: Short 
: Open 
: =50% 
: -20% 
: +20% 


>: Open 
: Open 


: +20% 


> +50% 
> =-20% 
: -202 
: +20% 
: +20% 
> +20% 
> +20% 
: Short 


OO Pwo 





hits or possible hits could have been accomplished. Thus, using the 
conventional procedure, within tolerance variations introduced 17% to 


42% errors. 


B. RESULTS USING SENSITIVITY FUNCTIONS 

Using the frequencies determined from the sensitivity functions a 
Set of signatures was generated for the six circuit conditions. In order 
to check the effect of variations within tolerance for this method the 
Same procedure as used in VI-A was followed. The result was 50 hits, 16 
possible hits and 6 misses. The 6 misses included 3 due to the test 
Signature showing normal response from the circuit. These results showed 
an error rate of 11% - 31%. 
Eight additional circuit conditions were simulated for each element 
to see what could be expected from a less artificial situation. Table 
VI-2. lists all the conditions set and the numerical result. The numerical 
result was assigned as follows: 


- closest library signature selected 

- right element selected 

- one of multiple signature closest to given condition 
-.one of multiple signatures right element 


mm & W Mm — 


- incorrect selection. 

ute, GISSEE and possible hits were also assigned as before. The 
result was 115 hits, 43 possible hits and 10 misses. Or, 68% hits and 
94% hits or possible hits. 


C. USING SUBROUTINE DET FOR CORRELATION 
___The same set of test signatures used in VI-B was used to test the 
programmed correlation and selection. The input and results are listed 


in Table VI-3. The results are coded as follows: 
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Run | Run 2 


Circuit 

Condition Vo/V, Vi/1, Vo/V, Vi/1, 

R1: Short ] 1 ] 1 
-~90% ] 1 ] 1 
~70% 1 1 ] 1 
-50% 1 1 1 ] 
~40% 1 1 1 1 
-20% 1 5 3 1 
+20% 1 3 1 3 
ai, 1 1 1 1 
+50% ] 2 ] 1 
+60% 1 1 ] 1 
+90% 1 ‘| ] 1 
x3 ] 1 1 1 
X75 1 1 } 1 
Open 1 1 ] 1 

l= short ] 1 ] 1 
-~90% 1 1 ] 1 
~70% || 1 2 2 
-50% 1 1 ] 1 
-~40% 1 1 ] 1 
-~20% 1 1 1 1 
+20% ] 5 ] 5 
+302 ] 3 ] ] 
+507 1 1 1 1 
+60% 1 1 ] 1 
+90% 1 1 ] 1 
X3 1 1 1 1 
75 1 1 } 1 
Open 1 1 1 1 

R2: Short 3 3 3 3 
~90% 3 5 3 2 
-/0% 3 5 3 5 
~50% 1 1 1 1 
~40% 1 1 1 1 
-20% 1 1 ] 1 
+20% 1 1 ] 1 
+307, 1 1 ie 1 
+50% 1 1 ] 1 
+60% 1 1 1 1 
+907 1 1 ] ] 
Open 1 ] ] 1 


TABLE VI-2. Results Using Sensitivity Functions. 
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correct primary selection 
one of two primary selections count 


WwW MP — 
] 


one of two secondary (no primary selections) 
4 - incorrect primary selection. 


The results were 61 rated 1, 12 rated 2, 21 rated 3 and 2 rated 4. 
Or, 63% completely correct and 98% with one of two selections correct. 
This is roughly the same result using manual correlation. The algorithm 
used for determining primary and secondary selection is outlined in 


Section VII-C. 
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VII. COMPUTER PROGRAMS 


Several computer programs were written in Fortran IV to accomplish 
much of the mathematical work associated with this study. Four of the 
programs are included herein with detailed descriptions. Each of these 
programs were used with several modifications. However, only the basic 
program is presented in each instance. All listings are in the designated 


‘section. 


A. PRIMARY ANALYSIS PROGRAM 

The primary analysis program contains the main program and three 
Subroutines. A fourth subroutine, PRQD, is used for root finding and 
1s not included in the listing. 

The program calculates the magnitude of both network functions as a 
function of the complex frequency, s, finds the roots of the polynomials 
involved and selects test frequencies. 

Each element in turn is then set to each out-of-tolerance value as 
selected previously, the magnitude of each function at each test fre- 
quency is calculated and all responses are quantized into signatures. 

The output consists of both network functions, the test frequencies, 
the nominal response at each test frequency and the list of library 
Signatures. 

The listing is the program which calculates the test frequencies as 
explained in Section IV-A. To modify the program for the selection from 
sensitivity functions the portion of the program between indicated com- 
ment cards was replaced with a read statement and the test frequencies 
were used as the input. 


Figure 21 is the flow graph for the main analysis program. 
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Figure 21. Flowgraph of Primary Analysis Program 
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B. TWO COMPONENT VARIATION PROGRAM 

The main analysis program was modified to generate signatures for 
Simultaneous two-component failure. Two different modifications were 
made, one generated signatures for only short and open element failures. 
The other allowed for short, -50%, +50% and open failures. 

Since the programs were only modifications to the main analysis 


program flow data are not included and the listings are omitted. 


C. FAULTY COMPONENT SELECTION 

This program was written to determine possible and probable component 
failures from the two input signatures. The signature from Vi/1, and 
Vo/N, are both processed. A component is flagged as possibly failed if 
one of the library signatures for that component is selected as a match. 
If two or more library signatures for a component are selected as a match 
that component is flagged as probably failed. 

The flow graph for the program is shown in Figures 22-24. The main 
program allows for several signatures to be read tn and sent to subroutine 
DET, one pair at a time. 

Subroutine DET first checks to see if the signature indicates the 
circuit is functioning properly (i.e., signature is all 5s). If so, a 
message is printed and control is returned to the main program. The first 
entry into DET with a non-nominal signature results in the library signa- 
tures being read in and stored. 

When it has been determined that non-nominal signatures are being 
processed, both input signatures are compared with the library signatures 
to determine the closest match distance. The signatures are compared a 


second time and those library signatures that are the minimum distance 
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Figure 22. 
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Figure 24. Flowgraph of Subroutine DET (Continued) 





from the input are used to set flags for determining possible and 
probable faulty components as described earlier. 
The program is listed showing the nearest neighbor matching 


procedure. 


D. SUBROUTINE LSQ 

Subroutine LSQ was written to determine which library signatures 
most closely matched the input signature. LSQ is similar in instruction 
flow to subroutine DET and therefore a flow chart is not presented. 

In order to determine a match, LSQ uses a method similar to the 
nearest neighbor concept. However, tf one considers the signature to 
be a vector with each test frequency a dimension, LSQ does not weigh 
every dimension equally. If one dimension can be strongly affected by 
random variations in component values that dimension is not considered 
equally with the other dimensions. Thus, the distance between the input 


vector and the library vectors becomes: 


ims 


2 
a. (X, - Y,) 
where 


X = bX, + boX5 + baX_ + ++: 


input vector 


Y = cyY, “ Ca¥, + c3Y, Be Paris 


library vector 


a. = weighting factor for the it? dimension (or 
frequency ) 
a. < |. 





LSQ operates on only one signature at a time and selects the proper 
library signature to match the input. The output is a list of the Te 
sible signatures and the component variation defined by each one. 

One may choose to set each le les This would be the case if the 
quantization process was sufficient to take care of the random variation 


problem. 
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VIII. SUGGESTED FURTHER STUDIES 


A. ALTERNATE APPROACH TO SIGNATURE TABLE 

Table VIII-1 lists the calculated response at each test frequency 
for the primary test circuit. Table VIII-2 lists the range of the 
response at each frequency and the maximum deviation from nominal response. 


Zero output is not included. 


Test Frequency (Hz) Response 
Vo/Vy Vi/T, 

00048 -2.28 Db 12.86 Ohms 

BZ, -6.77 Db 3.22 Ohms 

5055 -10.85 Db 14.32 Ohms 

299 -56.8 Db 3.01 Ohms 

2.99 -120.3 Db 4.34 Ohms 
49.97 =219.0° UD 53.46 Ohms 


TABLE VIII-1. Nominal Response at Test Frequencies 


Test Frequency (Hz) 


Response Range Maximum Deviation 


V/V, (db) Vj/1, (2)  Vo/Vy (db) Vy /1, (2) 

00048 0028 -29.91 3.0 66.77 -27.63 +53.9] 
.027 15.58 -46.16 .246 300.2 -39.39 +296,98 
055 -4 36 -55.7 0 > Soles ~44.85  +287.68 
299 -20.94 -97.14 p23 COG: ~40.34  +296.99 
2.99 -52.28 -160.41 3. 319.6  +68.02  +315.26 
49.97 -98.48  -255.5 3. 5338. +117.02  +5284.54 


TABLE VIII-2. Response Range and Maximum Deviation 


One can see that the responses do not cover the same range of values nor 
are the maximum deviation from the nominal response the same. From the nature 
of the network functions it can be assumed that the responses at. different 
test frequencies are not independent for a given circuit condition. Due to 


these situations it would seem logical to weigh equal changes at. different 
frequencies differently. 
5 ae oon es 50 





The Setting of quantization levels and subsequent signature 
generation is one way to accomplish this weighing. The setting of 
quantization levels based on the selected situations tends to weigh 
responses at some test frequencies less than others. However,. this 
is purely an "educated guess" method and has no analytical foundation. 

An alternate scheme was studied briefly and is suggested as the 
subject of further study. 

If one considers each test frequency as one of the dimensions in 
an n-dimensional space then each response would be a point in the n- 
dimensionai space. ineread of limiting calculations to a few selected 
Circuit conditions one could generate many points representing many 
circuit conditions. One could then have a large enough sample of data 
to do some statistical operations. 

In order to remove dependence the normal procedure is to multiply 
each vector, representing a point, by the inverse of the co-variance 
matrix. The result is the collection of points in a "transformed space" 
where dependency is accounted for. 

The problem is the construction of the co-variance matrix.. However, 
Since the data base is now sufficiently large,an estimate of the matrix 
can be calculated from the data. 

With a circuit under test the response would constitute a vector 
which would be multiplied by the inverse of the covariance matrix. At 
this point there are several options for fault identification. One can 
make a comparison with the stored data and make the fault identification 
using the nearest neighbor technique. An alternate approach would be to 
use an average of the nearest m neighbors when m is some positive integer. 


There are other possible procedures for matching. 
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If one is only interested in which component is faulty it might 
also be of interest to investigate any possible grouping of points which 
would be useful as in pattern recognition. The unknown could then be 
assigned to a group according to its position. Reference 5 details the 
necessary pattern recognition information. 

The above described procedure would require a computer with sufficient 
storage to process and store the information pertaining to the data base 
and covariance matrix. However such a computer is readily available for 


a nominal cost. 


B. EXTENSION OF PROCEDURE 

It soon becomes obvious that for even moderate sized circuits the 
procedure would be very difficult to use. For example, a 25-element 
circuit with 6 possible element value perturbations would yield 150 
Signatures. For ease in testing the number of test frequencies should 
be as small as possible. However, with 150 signatures considerable 
signature duplication is probable unless the number of test frequencies 
is large. 

Reference 3 shows how a ladder network may be divided into subnet- 
works such that the voltage gain of a given subnetwork is independent of 
the voltage gain of all preceding networks. Therefore, testing by add- 
ing subnetworks one at a time can yield the fault component. 

This procedure uses only the voltage gain but does show one approach 
to handing large circuits. In any case, a large network would have to 
be divided into subnetworks by some scheme in order to reduce the problem 


to a manageable size. 
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IX. CONCLUSIONS 


The majority of the conclusions made as a result of this study were 
based on a comparison of the results is pation VI and results published 
by other investigators. 

Most investigators reported approximately 75% isolation using one 
network function with no allowance for variation-within-tolerance of the 
-nonfaulty elements. One investigator [2] reported 78%-98% success. How- 
ever, the results were not sufficiently documented to judge the merit of 

the claim. | 

The results of this study show that an additional error rate of 
17%-42% can be introduced by the within tolerance variations if the 
conventional procedures are followed. However, if sensitivity functions 
are used to aid in selecting test frequencies the error rate due to other 
element variations drops to 11%-31% for one network function. Further- 
more, using two network functions fault isolation is successful for up 
to 94% of random single fault conditions. 

The high degree of successful fault isolation clearly points out the 
advisability of using at least two network functions for testing. The 
two functions used for this study were Vo/V) and V,/1,. However, the 
designed purpose of the circuit may dictate a different choice of 
functions to be used. . 

The results clearly show the improved accuracy of isolation using 
Sensitivity functions as a guide in test frequency selection. In addition 
to improved isolation this method also eliminates the possibility of 


unknowingly picking test frequencies which yield the same information. 
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